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THE  ATMOSPHERE  POLLUTION  CAUSED  BY  THE  TECHNOLOGICAL  LIQUID 
•GASES  EVAPORATION 

Zbigniew  Hulewicz,  Zygmunt  Jacklewski  and  Kamil  Eksanow 

The  calculation  method  of  the  atmospheric  dispersion  of  the  pollution, 
caused  by  the  liquid  gas  evaporation,  when  this  process  takes  place  in 
the  open  space,  has  been  discussed  and  illustrated. 

1 .  INTRODUCTION 

The  problem  of  the  gas  dispersion  in  the  atmosphere  during  the  present 
intensive  development  of  the  world  economy  becomes  a  very  important 
problem  involving  consideration  of  numerous  specialists  from  different 
branches  of  sciences  and  life.  This  problem  comprises  also  systematical 
and  willful  gas  ejection  into  the  atmosphere  as  well  as  unpredictable 
gas  penetration  into  the  atmosphere  as  the  result  of  reservoir  damages, 
industrial  installations,  etc. 

The  problem  of  systematical  and  controlled  gas  exhaust  into  the 
atmosphere,  mostly  from  industrial  chimneys,  is  theoretically  pretty 
well  worked  out  giving  details  connected  with  this  problem.  The 
computational  methods  of  this  type  of  the  pollution  spread  is  known. 
However,  the  scientists  are  short  on  the  descriptions  concerning  the 
cases  of  emergency,  when  the  liquid  gases,  especially  the  toxic  ones, 
which  dispersing  in  the  atmosphere,  can  cause  the  life  threat  even  in 
the  places  located  far  away  from  the  source  of  dispersion. 

The  present  study  gives  the  sketch  describing  the  method  of  the 
computation  of  the  gaseous  cloud  dimensions  as  well  as  the  area, 
where  the  life  threat  can  exist.  This  method  can  be  useful  for 
drafting  the  threat  maps  in  the  area  of  the  existing  installations 
or  the  installations  under  construction  or  in  the  area  with  the  liquid 
gas  reservoirs  depending  on  the  meteorological,  and  regional  conditions 
and  on  the  sort  of  gas.  This  method  can  serve  for  the  evaluation  of 
the  precise  location  of  new  constructions. 

2.  FUNDAMENTAL  PRINCIPLES 

The  case  of  the  reservoir  damage  with  the  liquid  gas  is  being  considered. 


Thi'-.  damage  caused  the  whole  amount  of  the  liquid  gas  being  ejected  outside 
However,  the  ejection  time  was  so  short  that  it  can  be  neglected.  It  is 
presumed  that  the  liquid  gas,  thrown  out  of  the  reservoir,  makes  "a  blot" 
having  round  shape  and  coating  uniform  thickness  all  over  the  surface 
and  the  surface  coated  with  the  liquid  gas  is  horizontal  and  non-absorbable. 
The  liquid  gas  blot  which  was  formed  in  this  way  can  be  considered  as  an 
intermittently  continuous  pointed  source  of  emission.  The  source  output 
depends  on  the  evaporation  speed  of  the  liquid  gar:  which  is  presumed  to  be 
constant. 

The  gaseous  pollution  cloud,  resulting  from  the  liquid  gas  evaporation, 
undergoes  diffusion  in  the  atmosphere  according  to  the  vorticity  diffusion 
laws.  After  the  liquid  gas  is  evaporated  completely,  the  gaseous  cloud 
keeps  moving  in  the  direction  of  the  wind  and  with  a  wind  speed  until  its 
complete  diffusion.  Cf  the  main  interest  is  the  area  where  the  gaseous 
concentrations  exceeding  the  permissible  ones  are  being  observed. 

The  problem  of  the  vorticity  diffusion  in  the  atmosphere  ground  layer  is  a 
complex  problem,  comprising  not  only  meteorological  factors  but  also 
tonographical.  Tr  order  tv  simplify  the  calculations  it  is  presumed  that 
during  the  process  of  the  liquid  gas  evaporation  and  spreading  the  gaseous 
cloud, the  meteorological  conditions  remain  unchanged  and  the  process  is 
taking  place  in  the  plane  area  which  is  partially  covered  with  trees  and 
partially  constructed. 

The  solution  of  the  problem  requires  considering  three  separate  processes, 
namely  -  the  liquid  gas  evaporation,  the  evaporated  gas  spread  at  the 
moment  of  evaporation,  that  is  from  the  emis.  ion  continuous  source  and 
the  gas  spread  when  the  evaporation  is  over. 

3.  THo  SP  IED  dF  THE  LIQUID  GAS  ORATION 

On  the  basis  of  the  boundary  layer  theory  Leibenson  /I/  has  introduced  the 
empirical  equation  for  the  liquid  evaporation  speed  from  the  free  surface  in 
the  open  space  which  ha  the  following  form: 
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where  q  -  liquid  evaporation  speed  /g/s/; 

D  -  diffusion  coefficient  of  the  liquid  steam  in  the  air  /cmVs/; 

M  -  molecular  mass  of  trie  liquid  /  g/mol  /; 

P  -  pressure  of  the  liquid  saturated  steam  under  the  Tq  temperature  /mm/Hg  /; 
T  -  surrounding  temperature  /K/; 

O  j  "Z. 

R  -  ga.  constant  which  is  equal  to  6. 237. cm  •  mmHg/K.mol; 
b  -  width  of  the  surface  of  the  evaoorating  liquid  /  cm  /; 

R  -  Keynolds  number  (R  =ul  / 0) ; 

G  6  P 

u  -  wind  speed  /  cm/s  /; 

1  -  length  of  the  surface  of  the  evaporating  liquid  /  cm  /; 

P  2 

j  -  kinematic  viscosity  of  the  liquid  stem  /  cm  /s  /; 

pr~  Prandtl  number  (Pr=\)/D);  0.265  -  experimental  coefficient. 

If  the  evaporating  liquid  surface  has  round  shape  ”'i til  the  d  diameter, 
then  the  following  dependence  is  compulsory: 
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Equation  /I/  is  not  being  solved  in  the  Reynolds  number  optional  scope. 
However,  it  is  possible  to  arrive  at  a  certain  generalization  of  equation  /I/, 
introducing  the  variable  quantities  of  the  k  experimental  coefficient  as  well 
as  these  of  the  n  rov.er  factor  under  the  Reynolds  number.  At  that  time 
this  equation  becomes: 
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Equation  /2/  can  be  transformed  to: 
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The  left  rid  of  thi~  equation  is  a  certain  dimen  ionless  quantity  designated 
as  A  and  in  this  case  one  obtains: 
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On  the  basis  of  equations.  /J/  and  /4/  it  becomes  possible  to  calculatef 
using  the  xp<.rimental  data  ,  the  values  of  the  A  dimensionless  number  for 
different  valuer  of  the  Reynolds  number.  The  A  -f(Re)  (tab.  1)  empirical 
dependence,  obtained  in  the  above-mentioned  way,  is  being  u'  ed  for  the 
calculation  of  the  liquid  evaporation  sreed: 
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Ci.  the  an -logy  of  the  End  law  of  the  Fick  diffusion,  the  changes  of 
the  gaseous  pollution  concentrate  -n  in  the  atmosphere  can  be  renresented 

by  the  following  equation: 


where : 
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e-  gar:  cone  on  rati  ; 

u , v , w  -  components  of  the  wind  speed  in  the  directi  in  of  the  !  ,Y,Z  axis 
of  the  co-ordinate  rectangular  system; 

k  ,k  ,k  -  com  onent  of  the  coefficient  of  the  vorticitv  diffusion  in  the 
x’  x’  z 

direction  of  th  X, Y,Z  axir  . 


The  princi  al  element  of  tee  the  <ry  of  th"  gaseous  pollution  dispersion  in 
th*.  tmo;  r  here  is  a  pro  es  interpretation  of  the  diffusion  coefficient  in  th 
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direction  of  the  Z  axis,  since  the  kind  of  surface, where  the  vorticity 

diffusion  takes  place,  influences  significantly  this  diffusion.  From  the 

[2] 

La jchtman^studies ,  it  turns  out  that  the  kind  of  the  basis  influences  first 
of  all  the  air  vortical  currents.  According  to  the  half empirical  theory  of 
vorticity,  the  kz  diffusion  coefficient  is  expressed  by  the  following 
equati  >n: 
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where  1  is  the  length  of  the  agitation  way. 


The  dependence  of  the  agitation  way  on  the  height  is  expressed  in  /3/ 


l  = 


<a> 


where;  y  -  constant  empirical  coefficien  (0,40); 

zq  -  parameter  characterizing  the  basis  roughness,  which  represents  the 
height  over  the  Earth  surface  where  the  wind  speed  is  equal  to  0  /m/ ; 
z  -  height  otoer  the  Earth  surface  /m/; 

(/  -  parameter  characterizing  the  degree  of  the  atmosphere  vertical  stability 
represented  by  the  following  function: 


In  our  geographical  latitude  the'b  values  fluctuate  from  -.3  uc  to  +.3, 
then: 


At 

e  <10> 

where: 

-  temperature  difference  at  the  50  and  200  cm  height  /°C/; 
u  -  wind  speed  at  the  z  height  /  m/s  /. 

Dependence  of  the  wind  seeed  on  the  height  /3,4/  is  represented  by  the 
following  equation: 

v  _  It i*/*»  (11) 
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where  u1  wind  speed  at  the  z1  height  /m/; 

After  the  transformation  one  can  obtain  from  equations  /7/,/8/,/11/  the 
following: 


The  expression  in  the  first  brackets  represents  the  k^  diffusion  coefficient 
in  the  vertical  direction  at  the  z^  height  and  in  this  case  one  can  write 


For  z.j=1  one  obtains: 

k,  -  M‘"‘.  (14) 

It  appears  from  this  that  in  order  to  calculate  the  coefficient, it  is 
necessary  to  know  the  k^  coefficient.  If  in  the  equation  determining  the 
k^  coefficient,  the  product  of  all  quantities  except  u^  is  designated  as  p, 
then  one  obtain 

(15) 

For  different  values  of  the  &  and  zQ  parameters  one  can  calculate  the  p 

values  and  from  them  it  is  possible  to  calculate  values  of  the  k^  coefficient. 

The  z  parameter  usage  is  complicated  because  it  changes  depending  on  £> ,  that 

is  on  the  atmosphere  condition.  Lajchtmn  /4/  has  introduced,  the  roughness 

parameter  de:  ignated  a'  z  which  depends  completely  on  the  basis  condition 

he  also  cslculatcd  the  values  of  the  /u^  quotient  for  different  vaiUo^ 

Fj  and  z  .  The  z  and  z  oerameter?  are  connected  between  each  other  with 
w  oo  o  00' 

the  following  de  endences: 

zo=zoo  in  the  c- re  of  isothermality ; 

zM=2z  in  the  case  of  convection; 

0  oo 

z  z  in  the  case  of  inversion 

o  oo 

Table  2  nresents  different  values  of  the  z  parameter  for  different  kinds  of 

oo 

the  basis.  Fig.  1  presents  the  nomogram  of  the  p  kxlut  -/(«,**.)•  funct'  -n. 
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Table  2.  Dependence  of  roughness  parameter  vs.  the 
kind  of  covering  of  the  terrain 


1 

Kind  ol'  covering  of 

! 

j 

Height  ! 

Roughness 

the  terrain  j 

(thickness) 

parameter 

1 

of  covering 

1 

1 

i 

l'"l  ! 

-  -  | 

[m] 

i 

t'losc-land -surface  with 

1 

1 

0,010 

no  plants 

j 

..  ..  _  _ _ _ 

0,05-0,1  1 

0,022 

Idle  land 

i 

0,3  j 

0,039 

i 

0,2 

0,032 

iil'O&i 

0,1  1 

0,025 

1 

0,05 

i 

0,015 

1 

1 

1 

1,2-1, 3 

0,046 

Itipo  corn 

0,060 

T.ipial  snow-covering  of 

1 

considerable  thickness 

! 

i 

0,010 

1 1  ll<'l|ll;ll  HtlUW-COVclillg 

of  average  tliiokuess 

i 

1 

'  do  10 

0,100 

Mingle  buildings,  clumps 

of  tiuslu  s  ami  true*. 

small  woods 

7 


The  vorticity  diffusion  coefficients  in  the  k  and  k  horizontal 

x  y 

surface  have  the  same  value  and  are  independent  from  the  height. 

They  are  designated  by  the  kQ  joint  symbol.  As  opposed  to  the 

vertical  direction,  the  diffusion  in  the  horizontal  direction  is 

not  a  subject  to  the  basis  hampering  influence.  Therefore,  at 

small  heights  the  k  coefficient  has  to  be  smaller  than  the  k 

z  o 

coefficient.  The  k  coefficient  is  increasing  along  with  the 

z 

height  increase,  whereas  the  kQ  coefficient  remains  constant.  At 
a  certain  height  in  the  layer  which  is  called  the  isotropic 
vorticity  layer  one  can  see  the  k z  and  kQ  coefficient  equalization. 
It  is  assumed  that  the  isotropic  vorticity  appears  at  the  height 
where  the  relation  of  the  squares  of  the  horizontal  component  of 
the  wind  speed  to  the  squares  of  the  vertical  one  reaches  the 
minimum  /5/.  The  xperience  indicated  that  for  the  isotherm,  the 
isotropic  vorticity  takes  place  at  the  12-14  m  height  /6/. 

Generally,  it  is  assumed  that  the  isotropic  vorticity  height  is 
equal  to  13m.  Using  this  assumption  one  can  calculate  the  kQ 
diffusion  coefficient,  which  is: 

kt  =  jfcj-is'-*.  (Hi) 

This  relation  can  be  used  in  the  inversion  condition 

for  a  small  wind  speed. 

Knowing  kz  and  kQ  one  can  solve  the  differential  equation  /6/ 
which  for  this  purpose  is  better  to  transform  introducing  kz 

according  to  formula  /13/  as  well  as  from  k  =  k  =  k  one  obtains: 

x  y  o 
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Equation  /7/  is  being  solved  while  maintaining  the  following  boundary 
conditions : 


a)  it  is  assumed  that  there  is  no  gas  absorption  through  the  basis,  that  is 


do 

lun  —  =  0 : 

e  OZ 


b)  at  the  t  -  -  0  time  the  gas  concentration  in  all  points  of  the  space  is 
equar  to  z». ro,  except  for  the  initial  noint  of  the  co-ordinate  system 
where  there  is  the  source: 

linw’(r,  X,  y,  z)  0  for  T  “  °»  0,  *  -  0; 

c)  at  the  distance  infinitely  remote  from  the  source,  the  gas  concentration 
is  always  equal  to  zero: 

liiuo(r,  a),  y,  *)  =  0,  r  --■=  f 

r  »oo 

the  to:  '  r  'a -it  •>  th  r  •  role-'  •  e«  through  trie  "urce  is  located  in  the 
atmosphere: 

v  »».-/*  'f  a, 

“  -00  -no 

where:  q  -  the  emission  sourse  output  is  equal  to  the  evaporation  speed  of 
the  liquid  gas,  :/g/s/; 

^  -  evaporation  time  of  the-  liquid  gas  /s/; 

c(x,x,  y,z)  —  f  '  c  -.centration  at  the  point  of  the  x,y,z  co-ordinates  after  -fki 
time  '  -  c '  f.  ornent  of  evaporation  beginning  /g/m  /. 

The  solution  of  equation  /I  7/  after  taking  into  consideration  all  above- 
mentioned  boundary  conditi'ns  will  look  like: 


'‘(J-.K.  I 
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where:  /«  1  |  *.  '  the  index  of  the  atmosphere  stability;  /’(l-pl/w) 

C  function  of  the  l  \  I  jm .  -gument; 
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Assuming  that  z  =1  ,  equation  /1 8/  will  be  reduced  to 

q(nilklm')'lm  r  u  jyt  Um  i  , 
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the  following: 
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Equation  /1 9/  enables  to  calculate  the  gas  concentration  noting  from 

the  emission  sourse  with  the  q  constant  output.  This  calculation  can  be 
made  at  any  roint  with  the  x,y,z  co-ordinates. 

4.  GAS  CLOUD  DIFFUSION  A FT EH  EVAPORATION  IS  OVER. 


The  gas  cloud,  created  as  the  result  of  the  liquid  gar  evaporation  after 
the  reservoir  or  installation  was  damaged,  will  detach  from  the  emission 
source  after  the  bp  time  which  ir  equal  to  the  evanoration  time.  After 
the  emission  is  over,  the  gar  cloud  will  keep  diffusing  as  the  result  of 
stormy  atmospheric  diffusion. 

The  formula  for  the  calculation  of  the  gar  in  atmospher  at  different 

\  /T' 

distances  from  the  emission  sourcepr  after  che  optional  V  time  expiration, 
can  be  introduced  by  means  of  the  integration  of  equation  /I  7/  while 
maintaining  the  a,  b,  c  boundary  conditions  as  well  as  the  approximative 
condition  which  presumes  that  the  total  amount  of  the  liquid  gar  after 
evanoration  com  roses  th  cloud  with  complete  dimensions. 

quation  /I?/  after  being  solved  ^assuming  that  z^l)  aquires  the  following 
form  in  this  care: 

x  ‘  I  y  -  v 

klmtTu  +  h',il 

1  \  ,  *  (20) 

where:  Q  -amount  of  gas  in  the  cloud  /g/;  ^u  -  time  measured  since  the 
moment  of  the  evaporation  end  /s/;  r,h  -  specific  dim- ns ions  of  the  gas 
cloud  (r-  radius  of  toe  horizontal  profile  on  the  Earth  level; 
h  -  height  /m/) ; 

The  gas  cloud  moves  in  the  direction  of  the  wind  with  its  speed.  In  the 
same  way  tie  co-ordinate  system  where  eouation  /20/  is  mandatory,  has  to 
move  in  the  a  me  -'-y.  It  is  assumed  that  the  beginning  of  the  co-ordinate 
system  is  cover'd  with  the  gravity  center  of  toe  gas  cloud. 

The  r  and  h  parameters  introduced  in  equation  /?0/  are  not  absolute 
dim-nr i  >n-  of  th-  gas  cloud.  Jn  order  to  determine  the  values  of  these 
naramrt  rs,  eouati''n  /?')/  can  be  1  ogarithmated  wich  gives  the  following 
resul t : 
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*/’  1  + 
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111  6-  —  lu  ^4  — 


x*  \  -y* 


where 


nr 


4fc«Tw  -f  r*  k  t  m*Tu  )-*“  ’ 

_ Q 

|l  |  -|-J  (A:,to*tu  +  (4A:0tm  1  r») 


If  the  cl oii:i  ■■  ecti  >n  i:  being  regarded  in  the  z=  const.  plane,  then 


.'here 


x*  I  i/‘ 


J‘  *,»4. 1-*-)’ 


(22) 


Equation  /Z?J  can  be  written  in  the  form: 

.  Ah 

xM-0*  -  (41.t„  f-r*)!!! - . 


(23) 


This  is  the  circle  equation  with  the  radius : 


(41«th  |  r*)ln 


AH 


(24) 


«  1  p  p 

From  equation  /?4/  it  appears  that  if  Ui=0  and  ABc~  =e,  then  r  =R  . 

_  1  _  p 

If  in  this  case  z=h,  then  B=e  and  in  the  result  e=Ae  Therefore,  r  and  h 
are  such  specific  cloud  dimensions  for  which  the  gas  concentration  in  the 
atmosphere  at  the  tu=0  moment  amounts  to: 


(J 


(2f») 


5.  GBIToRAL  METHOD  OF  CALCUL  IT  I  "NT, 

The  basic  data  for  calculations  to  be  determined  are  the  parameters 

1  1 


diameter izing  the  liquid  gas,  such  as:  the  ga?  sort  and  amount,  molecular 
mass,  pressure  of  the  gas  saturated  under  the  temperature  for  which 
calculations  will  be  done,  density,  kinematic  viscosity,  the  diffusion 
coefficient  under  the  same  temperature,  permissible  gas  concentration  in 
the  air;  the  parameters  characterizing  the  atmosphere  condition:  the  wind 
s  eed  and  direction,  the  air  temnerature  at  the  50  and  200  cm  height; 
and  finally  the  parameter  of  the  basis  roughness  designated  as  zQQ-  It 
is  necessary  to  evaluate  the  size  of  the  liquid  gas  blot. 

The  gas  physical  and  chemical  parameters  are  stabilized  and  given  for  the 
standard  temperatures.  Values  of  these  parameters  for  ;c:ier  temperatures 
can  be  calculated  from  the  formulas  given  in  th>"  lite  ature. 

After  the  determination  of  the  basic  data,  it  becomes  possible  to  calculate 
the  Reynolds  number,  as  well  as  the  Prandtl  number  and  A  number  from 
formula/5  /  which  is  a  formula  of  the  evaporation  speed  of  the  liquid  gas 
which  determines  at  the  same  time  the  output  of  the  emission  source. 

The  calculation  of  the  size  of  the  area  where  the  gas  concentration  in  the 

air  e  coeds  the  permissible  concentration  lies  in  the  determination  from 

equation  /1 9/  the  x,y,z  coefficients  of  the  space  points  where 

c ,  .  =c,  (c,  -  permissible  concents  tion).  The  y  co-ordinates  are  being 

vXjYjZ/  g  a 

calculated  for  the  following  x  co-ordinates  (e.g.  the  pitch  is  every  10  m) 

/fj 

within  the  scale  from  x=0  un  to  x=U  u  (  v  -  evaoorati  n  time;  u  -  wind 

p  p  ’ 

speed  at  a  given  height).  The  calculations  are  being  conducted  for 

differ  nt  z  heights  or  for  the  chosen  height.  However,  the  most  frequently 
chosen  height  is  2  meters,  since  in  this  atmospheric  layer  the  human  being 
is  threatened  .nos  tly. 

The  coordinates  of  the  points  where  c=c^»  calculated  in  this  way,  denote 
the  cloud  dimensions  at  the  moment  when  the  eva  oration  is  over.  After  the 
evaporation  is  over,  the  cloud  is  moving  in  the  direction  of  thr  "ind. 

The  cloud  location  and  dimensions  after  the  (/  optional  time,  calculated 
from  the  moment  of  the  evaporation  termination,  are  being  calculated  from 
f  rmula  /20/.  However,  there  are  two  problems,  whose  soluti  n,  is  conditioned 
by  a  correct  interpretation  of  the  calculation  results.  The  first  problem 
is  the  d  termination  of  the  r  and  h  specific  dimens  ions  of  the  gaseous  cloud, 
that  is  the  parameters  ■  hich  are  necessary  when  equation  /20/  is  being  used. 
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The  s  cond  problem  lies  in  the  non-congruence  of  the  figures  obtained  as 
the  result  of  the  equation  /1 9/  and  /20/  solution. 

It  is  proposed  to  use  the  following  procedure  which  represents  a  conditional, 
approximative  solution  of  the  above-mentioned  problems:  the  figure  vhich  is 
being  madd  by  the  c=const.  lines,  after  the  0  ^time  (solution  of  equation  /19/) 
should  be  tran  ferred  in  proper  dimension  of  the  circles.  The  transposition 
basis  can  be  maintaining  the  same  area  of  the  surface,  circum.-cfcibed  by  the 
given  figure  and  circle.  After  obtaining  two  circles  with  the  and  R? 
radius  which  correspond  to  .he  c1  and  c^  concentrations , it  becomes 
possible  to  determine  the  specific  radius  of  the  r-cloud  ueJ:.&  equation  /?4/. 
The  specific  height  h  can  be  calculated  in  the  ••me  way. 

<r 'J'sj 

In  order  to  calculate  the  dimensions  of  t‘e  cloud  after^Gutime  it  is 
possible  t  apply  equation  /'4/  or  /?  /.  The  circle  obt  ined,as  the  result 
of  calculations  ith  the  R  radius  has  to  be  retransferred  (maintaining 
the  same  size  of  the  surface  area)  into  the  figure  similar  to  the  riginal 
one. 

This  procedure  enable  to  determine  not  only  the  size  of  the  gaseous  cloud 
but  also  its  shape.  However,  one  must  be  warned  that  the  assumption  of  the 
similarity  of  the  cloud  share  at  the  moment  of  its  diffusion  not  always 
can  be  realized. 

6.  z.  AM.  LARY  calcut,'\tton:; 


The  mentioned  method  shows  the  calculations  of  the  results  of  the  damage  of 
the  reservoir  which  contains  the  liquid  chlorine.  The  calculations  were 
made  by  me  ns.  of  the  Minsk  22  and  Odra  1004  digital  comnuters.  The  basic 
data,  are  compared  in  Table  3  and  4. 

The  Re,  A  values  (Tab.1)  and  q  values  (Formula  /5/)  calculated  on  the  basis  of 
these  data,  armunt:Re=  8.0*  10^;  A=2.5*1o\  the  evaporation  speed 
q=37.7  kg/s. 

It  is  assumed  that  the  process  of  the  atmor-pheric  diffusion  takes  rlace  in 
the  inversion  conditions,  in  the  area  which  is  partially  covered  with  forests 
and  partially  cultivated.  For  the  inversion  the  average  value  of  the  ts 
parameter  ls»20  (m=1+6  =  1.2);  however,  the  surface  roughne:  s 


in  the 


accepted  here  area  amounts  .10m.  For  these  conditions  one  reads  out  the 

*  iUft 

value  of  the  k^/^  quotient  equal  to  .04  fj.om^monogram  (Fig.1).  From  this 
value  me  calculates  the  value  of  the  diffusion  coefficient  when 
assuming  that  the  wind  speed  at  the  ?m  height  is  \ .?  times  lover  than  this 
one  at  the  Pm  height,  that  is  u^=  u0/l .2.  These  data  are  comnared  in  Table  5. 

Using  formula  (I  .1)  the  chlorine  diffusion  in  the  atmosphere  has  been 
com-uted.  The  gaseous  cloud  dimensions  indicate  the  oints  "here  the 
chlorine  concentration  reaches  the  permissible  value.  It  is  assumed  that 
the  chlorine  permissible  concentration  amounts  .001mg/dm^.  'his  value  is 
given  by  the  Polish  State  Standard  PN-56/Z-040030  which  identifies  the 
maximum  rermir  ible  gas  concentrations  in  the  air  surrounding  the  v orking 
'lace.  The  gaseous  cloud  dimensions  after  the  va’ oration  v a  terminated  are 
presented  on  Fig.  ?. 


3.  Start  data  fur  computing 


Quantity 

Symbol 

1 

i 

!  Numerical 
value 

1 

Amount  of  chloride 

V 

1 

j  1000  kg 

S|>ol  (limur.trr 

d 

10  III 

l 

» 

Air  tiiopriatoro 

T 

203  K 

Wind  velocity  at  a  higlit  of  2mj 

“a 

4  m /a 

iVrn  Hum  Mo  ooiicoiitialion  ofOlg| 

i 

'A 

0,001  g/in* 

Toxic  concentration  of  cl,  j 

" 

0,012  g/nr* 

O.OiKf/m* 


Tublu  4.  1‘byoieo -chemical  properties  of  chlorine  iu  10°  C 


Quantity 

Symbol 

N  umeruul 
value 

Molecular  uiuiM 

M 

70, HI  ({/mol 

Saturated  vapour  procure 

4903  uim  lig 

Kino  tic  viMco^ity 

V 

0,0442  cm*/* 

Uifliuiuii  coefficient 

b 

0,124  cm*/* 

-  Table  S 


Quantity  j 

Symbol 

N  umerieal 
value 

’  *uo 

0,10  m 

Kougliuoxa  parameter 

Vortical  stability  |iaramoter 
of  tbo  atmosphere 

1 

!  * 

i 

| 

\  0,20 

Index  at  ability  of  tbo 
atmosphere 

m  -  1  f  ( 

1,20 

Diffuaiou  ooeffieieut  (eq.  (IS)) 

0,133  in*/a 

fc» 

1,035  m */• 

lliffuaion  coefficient  (eq.  (10)) 

i 
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Fig.  2a.  Vertical -hoc  Uou  XZ  and  cross-section  X  X  of  a  chlorine  cloud  for  the  inversion  at 

the  end  liuic  of  source  emission 


The  cloud  recti  on  at  the  moment  of  the  ev-an  oration  termination  in  the  z=?  m 

plane  har  the  shape  cl  ore  to  the  isosceles  triangle.  Two  triangles  were 

found.  On  ti.oir  sides  the--e  are  the  known  c^  and  Cp  concentration.  The  areas 

of  these  triangles'  surfaces  have  be' n  calculated,  -\fte  that  the  radius 

of  the  circles  whose  area  is  equal  to  the  triangles'  areas  have  been  found. 

For  calculation-  the  following  concentrations  have  been  accepted: 

c.=  .001  mg/dm^;  c_=  .?■  0  mg/dm^.  The  obtained  radiuses  were:  R, -43.50  m; 

1  1 

R.^42»  1  -m.  Using  formula  /?  /  as  a  basis  and  as*  using  that  l/yl  =0,  it  was 
obtained  a  specific  radius  of  the  cloud  -  r=1~.75ra. 

The  gareaur  cloud  recti  on  in  the  y-0  vertical  section  has  also  (it  is 
po  sible  to  assume  in  amroximation)  the  triangle  shaoe.  T  o  trianeles 
corresponding  the  c,and  c0-  concentrations  have  been  found.  Their  surface 
areas  ha v  been  d-lerminea  -aid  then  one  obtained  the  circle  sections 
having  the  same  areas  and  chard"  being  equal  to  the  RR^and  2P.  0  diameters. 

The  calculated  arrows  of  these  "ections  amount  to  z  =?7. oOm  and  z?-  14.55  m. 


Applying  these  data  it  has  been  obtained  from  equation  /--I/  with  the 

/r 

condition  that  l/=0  that  the  h-m'-ci  f ic  hei  g  .t  is  equal  to  3.0  3m. 

Knowing^r  and  h  vhic.a  ar*.  specific  parameters  of  th  cloud,  one  calculated  on 

the  basis  of  eouati  n  ( ?_U )  the  radiuses  and  the  following  areas  of  the 

u  =  2  ai  # 

j;o  ’izontal  recti  ns  of  tie  cloud  in  the.  -one  after  different  \J  time 

*  u 

for  two  concentration  values  of  th-  chlorine  in  tso  air,  that  is  for  the 


permissible  concentration  (.001  mg/d m')  and  the  concent- 


t:L 


c  n  l  s 


beyond  '-nduranc  of  the  human  being  even  for  short- lasti ng  expocare. 

Thi:  concentrati  >n  i-  evaluated  as  .01?  mg/dm'  /7/  and  called  "toxic"  and 
designated  a:  c^.  Th.  next  step  was  to  find  out  the  triangles  with  the  areas 
being  equal  to  the  area  of  tr-.o  horizontal  s- cti -ns  of  the  cloud.  These 
ttiangles  ire  similar  to  the  triangular  appearing  at  fhe^  time.  These 
triangles  represent  approx  im  :  t  ive  dimens  -nr  and  rhaso  of  the  chi  'trine 
cloud  at  the  in  ment  of  its  di  f  fusion.  Tney  are  sho  n  o-  Fig.  3.  It  is 
natural  that  the  gaseous  cloud  hiving  the  shape  a  pr uximating  the 
triangle  doer  not  possess  its  sharply  outlined  elemen.tr,  especially 
vertexes.  It  is  evamrlified  on  th--  nicture  by  mean-  of  likening  the  shape 
of  the  diffusing  cloud  to  the  z-?m  faction  on  Fig.?. 


The  change  of  t  c  cl  >ud  dimensi  n  at  to-  moment  of  its  dif-'u-ion  when  the 
emi  ion  is  -ver  i  r-'-rr- srnted  in  a  continuour  wai  on  Mg  4  as  the 
change  of  its  fr~n  -idth. 


!ne  erS6ntial  Valuei:  cabling  to  conduct  the  evaluation  of  th  tv, 
in  a  specified  distance  fro*  tv  •  evaiuotx  ,n  of  the  threat 

anCe  frcra  the  esus  ion  source  are  the  fnilA  • 
the  time  afte  -hi ch  th *  „  tile  follov:ing  ones: 

.■ni.cn  the  cloud  front  will  reach  thi-  _ 

during  "hi ch  this  ni  ce  -m  r-  •  pl  C0;  the  tiroe» 

-  .Ce  ’..-ill  remain  within  the  cloud  * 

times  on  the  distance  from  th*  ■  "  ‘  Th  dc-~endenc"  of  there 

fr0m  the  eolSKion  source  in  shown  on  Fig.5. 


Fig.  3.  Zone*  of  chioritia  rink 


0) 


Fig.  ».  Start -limo  ami  dinalmn  «f  i-.hWim>  risk  in  tin*  wind  direction  for  a  given  distance 

from  source  emission 


"»  -=  •-  "a  ■*  «.«U  f  «/iu*.  C|  -■  o.rtl-’  g/m3.  <J  ~=  1 000  kg,  U  =»  10  m,  u,  -  4  m/s,  4  -  2  m 


7.  DTFCU.:  ION 


The  method  Presented  hereof  calculation  of  the  ga.  diffusion  in  the  n.ear- 
Karth  layer  of  t:,o  atrao.onhe  e  in  not  a  perfect  method.  It  has  certain  faults 
resulting  from  the  acceptance  of  numerous  sin  lifying  as-umr ti-nn .  Such 
assumptions  are  the  statements  that  in  the  ca  e  of  damage  l e.g.  a  reservoir) 
a  liquid  gas  after  pouring  out  compose  a  blot  with  permanent  thickness; 
that  the  surface  of  tec  liquid  level  during  th  -.•va-orati  n  r  min.  remanent 
and  that  th*  eva  oration  -p-  ed  is  a  reman  nt  value,  ~n  ro-lity,  the 
eva  oration  so end  can  undergo  ignlficant  fluctuations  cau  ed  not  only  by 
a  possible  change  of  the  evaporated  s,urfa.c-  but  al  -a  by  the  bn  i  nb  -orption, 
by  the  wind  speed  fluctuation  ,  by  changes  of  th  t: ■  rnal  conditions. 

The  assumoti  sn  of  th  -  const -noy  of  th  -  evaoorati  >n  spued  vas ,  ho-  .-v-r,  toe 
condition  of  tri-  oluti  n  of  equati  n  /I  1/. 

The  mathematical  mo-  1  v/hich  is  represented  by  formula  /1 9/  is  also  idealized 
It  does  not  t-  e  into  consideration  such  f-ctors  as  th-'  r  cul  tur  of  the 
.varth 1 '  urfac-- ,  vro:  nc  -  of  dlf'e-ent  regi  >nul  h  :.d  ances,  changes  nf  the 


2"l 


atmospheric  condition."  a  well,  as  the  fluct  lation.*;  of  the  wind  speed  and 
direction. 

The  ram  remarks  are  related  to  the  calculation."  of  the  gn  eour  cloud 
diffu; ion  after  the  ova  oration  termination  (formula  /20/) .  When  introducing 
this  formula  it  ha"  been  a.  .  umed  that  th*  total  gas  amount  after  evaporation 
is  in  the  atmosphere,  whereas  a  part  of  the  gar  is  being  absorbed  by  the 
oil,  water  reservoir",  forests,  etc.  It  has  also  been  ar  sumed  that  there  is 
imilarity  of  t  o  go  eour,  cloud  rhaoe  during  the  evaporation,  which  not 
alvay  can  bo  realized. 

.tome  f  the  mentioned  simplifying  a  um-'tionr  worsen  the  results  of 
calculations,  others  act  in  a  different  directi  n.  it  is  difficult  to 
determine  vhich  ones  of  t  sea  will  prevail,  ao-ever,  one  can  assume  that 
these  action.'  can  be  com  enra'-ed  in  certain  conditions.  There:  ore,  it  is 
also  considered,  that  thi"  ra  >thod  ca  be  applied  to  th  purpose  for  which 
it  was  adopted,  that  is  for  the  calculatio  .s  of  the  result"  of  damage  of 
the  rerervoirr  • ith  liquid  gases  and  analogous  phenomena. 
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